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I. INTRODUCTION 

Tunneling is a supreme quantum effect. Every introductory text[T] on quantum mechanics gives the paradigm 
example of a particle tunneling through a one-dimensional potential barrier despite having a total energy less than 
the barrier height. Indeed, the reader typically works through a number of excercises, all involving one-dimensional 
potential barriers of one form or another modelling several key physical phenomena ranging from atom transfer re- 
actions to the decay of alpha particles [2] • However, one seldom encounters coupled multidimcnisonal tunneling in 
such texts since an analytical solution of the Schrodinger equation in such cases is not possible. Interestingly, the 
richness and complexity of the tunneling phenomenon manifest themselves in full glory in the case of multidimen- 
sional systemsfJJ. Thus, for instance, the usual one-dimensional expectation of increasing tunneling splittings as one 
approaches the barrier top from below is not necessarily true as soon as one couples another bound degree of freedom 
to the tunneling coordinate. In the context of molecular reaction dynamics, multidimensional tunneling can result in 
strong mode-specificity and fluctuations in the reaction rates |4j. In fact, a proper description of tunneling of electrons 
and hydrogen atoms is absolutely essential [SJ [3] even in molecular systems as large as enzymes and proteins. Although 
one usually assumes tunneling effects to be significant in molecules involving light atom transfers it is worth pointing 
out that neglecting the tunneling of even a heavy atom like carbon is the difference between a reaction occuring or 
not occuring. In particular, one can underestimate rates by nearly hundred orders of magnitude [7]. Interestingly, 
and perhaps paradoxically, several penetrating insights into the nature and mechanism of multidimensional barrier 
tunneling have been obtained from a phase space perspective[8l[9]. The contributions by Creagh, Shudo and Ikeda, 
and Takahashi in the present volume provide a detailed account of the latest advances in the phase space based 
understanding of multidimensional barrier tunneling. 

What happens if there are no coordinate space barriers? In other words, in situations wherein there are no static 
energetic barriers separating "reactants" from the "products" does one still have to be concerned about quantum 
tunneling? One such model potential is shown in Fig. [l] which will be discussed in the next section. Here we have the 
notion of reactants and products in a very general sense. So, for instance, in the context of a conformational reaction 
they might correspond to the several near-degenerate conformations of a specific molecule. Naively one might expect 
that tunneling has no consequences in such cases. However, studies over last several decades [TUHID] have revealed 
that things are not so straightforward. Despite the lack of static barriers, the dynamics of the system can generate 
barriers and quantum tunneling can occur through such dynamical barriers[2I]. This, of course, immediately implies 
that dynamical tunneling is a very rich and subtle phenomenon since the nature and number of barriers can vary 
appreciably with changes in the nature of the dynamics over the timescales of interest. This would also seem to 
imply that deciphering the mechanism of dynamical tunneling is a hopeless task as opposed to the static potential 
barrier case wherein elegant approximations to the tunneling rate and splittings can be written down. However, 
recent studies have clearly established that even in the case of dynamical tunneling it is possible to obtain very 
accurate approximations to the splittings and rate. In particular, it is now clear that unambiguous identification of 
the local dynamical barriers is possible only by a detailed study of the structure of the underlying classical phase 
space. The general picture that has emerged is that dynamical tunneling connects two or more classically disconnected 
regions in the phase space. More importantly, and perhaps ironically, the dynamical tunneling splittings and rates 
are extremely sensitive to the various phase space structures like nonlinear resonances [22 -29 , chaos [16-rT9l [30l [3T| 
and partial barriers [32l [33]. It is crucial to note that although purely quantum approaches can be formulated for 
obtaining the tunneling splittings, any mechanistic understanding requires a detailed understanding of the phase 
space topology. In this sense, the phenomenon of dynamical tunneling gets intimately linked to issues related to phase 
space transport. Thus, one now has the concept of resonance-assisted tunneling (RAT) and chaos-assisted tunneling 
(CAT) and realistic systems typically involve both the mechanisms. 

Since the appearance of the first book[5T] on the topic of interest more than a decade ago, there have been 
several beautiful experimental studies [MH5^ that have revealed various aspects of the phenomenon of dynamical 
tunneling. The most recent one by Chaudhury et al. realizes[lD] the paradigmatic kicked top model using cold ^'^■^Cs 
atoms and clearly demonstrate the dynamical tunneling occuring in the underlying phase space. Interestingly, good 
correspondence between the quantum dynamics and classical phase space structures is found despite the system being 
in a deep quantum regime. As another example, I mention the experimental observation 41^ by Foiling et al. of 
second order co-tunneling of interacting ultracold rubidium atoms in a double well trap. The similarities between 
this system and the studies on dynamical tunneling using molecular effective Hamiltonians is striking. In particular, 
the description of the cold atom study in terms of superexchange (qualitative and quantitative) is reminiscent of the 
early work by Stuchebrukhov and Marcus [12] on understanding the role of dynamical tunneling in the phenomenon 
of intramolecular energy flow. Further details on the experimental realizations can be found in this volume in the 
articles by Steck and Raizen, and Hensinger. An earlier review [43| provides extensive references to the experimental 
manifestations of dynamical tunneling in molecular systems in terms of spectroscopic signatures. Undoubtedly, in 
the coming years, one can expect several other experimental studies which will lead to a deeper understanding of 



dynamical tunneling and raise many intriguing issues related to the subject of classical-quantum correspondence. 

As remarked earlier, it seems ironic that a pure quantum effect like tunneling should bear the marks of the underlying 
classical phase space structures. However, it is useful to to recall the statement by Heller that tunneling is only 
meaningful with classical dynamics as the baseline. Thus, insights into the nature of the classical dynamics translates 
into a deeper mechanistic insight into the corresponding quantum dynamics. Indeed, one way of thinking about 
classical-quantum correspondence is that classical mechanics is providing us with the best possible "basis" to describe 
the quantum evolution. The wide range of contributions in this volume are a testimony to the richness of the 
phenomenon of dynamical tunneling and the utility of such a classical-quantum correspondence perspective. In this 
article I focus on the specific field of quantum control and show as to how dynamical tunneling can lead to useful 
insights into the control mechanism j44l |45]. The hope is that more such studies will eventually result in control 
strategies which are firmly rooted in the intutive classical world, yet accounting for the classically forbidden pathways 
and mechanisms in a consistent fashion. This is a tall order, and some may even argue as an unnecessary effort in 
these days of fast computers and smart and efficient algorithms to solve fairly high dimensional quantum dynamics. 
However, in this context, it is useful to remember the following which was written by Born, Heisenberg, and Jordan 
nearly eighty years ago '46^ 

" The starting point of our theoretical approach was the conviction that the difficulties that have been encountered at 
every step in quantum theory in the last few years could be surmounted only by establishing a mathematical system for 
the mechanics of atomic and electronic motions, which would have a unity and simplicity comparable with the system, 
of classical mechanics... further developement of the theory, an important task will lie in the closer investigation of the 
nature of this correspondence and in the description of the manner in which symbolic quantum geometry goes over 
into visualizable classical geometry." 

The above remark was made in an era when computers were nonexistent. Nevertheless, it is remarkably prescient 
since even in the present era one realizes the sheer difficulty in implementing an all-quantum dynamical study on 
even relatively small molecules [27|. In any case, it is not entirely unreasonable to argue that large scale quantum 
dynamical studies will still require some form of an implicit classical-quantum correspondence approach to grasp the 
underlying mechanistic details. With the above remark in mind I start things off by revisiting the original paper[TT] 
by Davis and Heller since, in my opinion, it is ideal from the pedagogical point of view. 

II. DAVIS-HELLER SYSTEM REVISITED 

Three decades ago, Davis and Heller in their pioneering studvfTT] gave a clear example of dynamical tunneling. 
A short recount of this work including the famous plots of the classical trajectories and the associated quantum 
eigenstates can be found in Heller's article in this volume. However, I revisit this model here in order to bring forth 
a couple of important points that seem to have been overlooked in subsequent works. First, the existence of another 
class of eigenstate pairs, called as circulating states[TT], which can exert considerable influence on the usual tunneling 
doublets at higher energies. Second, a remark in the original paperjTT] which can be considered as a harbinger for 
chaos-assisted tunneling. As shown below, there are features in the original model that are worth studying in some 
detail even after three decades since the original paper was published. 

The Hamiltonian of choice is the two degrees of freedom (2DoF, in what follows the acronymn DoF stands for 
degrees of freedom) Barbanis-like modeljll] 

H{s,u,Ps,Pu) - 2 (P? +P«) + 2 {^y + -^W) + ^su" (1) 

with the labels 's' and 'u' denoting the symmetric and unsymmetric stretch modes respectively. The above Hamiltonian 
has also been studied [JS] in great detail to uncover the correspondence between classical stability of the motion and 
quantum spectral features, wavef unctions, and energy transfer. The potential is symmetric with respect to u O —u 
as shown in Fig. [l] but there is no potential barrier. Davis and Heller used the parameter values ujg — 1-0, a;„ — 1.1, 
and A = —0.11 for which the dissociation energy E^is = w^w^/SA^ = 15.125. Note that the masses are taken to be 
unity and one is working in units such that h = \. The key observation by Davis and Heller was that despite the 
lack of any potential barriers several bound eigenstates came in symmetric-antisymmetric pairs \ipi) and IV'2) with 
energy splittings much smaller than the fundamental frequencies i.e., AE = \Ei — £2] -^ 0{1) . In Fig. [Tithe various 
splittings between adjacent eigenstates are shown and it is clear that several "tunneling" pairs appear above a certain 
threshold energy. 

How can one understand the onset of such near degeneracies in the system? The crucial insight that Davis and 
Heller provided was that the appearance of such doublets is correlated with the large scale changes in the classical 
phase space. The nature of the phase space with increasing total energy is shown in Fig. fusing the (u,p„) Poincare 
surface of section. Such a surface of section, following standard methods, is constructed by recording the points, with 
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FIG. 1. The left plot shows the two-dimensional potential used by Davis and Heller to illustrate the phenomenon of dynamical 
tunneling. Note that there are no barriers in the coordinate space. On the right the magnitude of the splitting between adjacent 
eigenstates are shown as a function of the eigenenergies. The first and second dotted vertical lines indicate the energies around 
which the symmetric (s) and unsymmetric (u) stretch modes become unstable. Sequences of doublets with very small splittings 
can be seen around the onset of the instabilities. On the other hand, the regularity of the sequences vanishes at higher energies. 



niomentum p^ > 0, of the intersection of a trajectory with the the s = plane in the phase space. Such a procedure 
for several trajectories with specific total energy generates a typical surface of section as shown in Fig. [2] and clearly 
indicates the global nature of the phase space. It is clear from the figure that the phase space for i5 < 11 is mostly 
regular while higher energy phase spaces exhibit mixed regular-chaotic dynamics. One of the most prominent change 
happens when the symmetric stretch periodic orbit {u = 0,p„ = 0) becomes unstable around E w 6.8. In Fig. ^the 
consequence of such a bifurcation can be clearly seen for .E = 9 as the formation of two classically disjoint regular 
regions. In fact, the two regular regions signal a 1:1 resonance between the stretching modes. The crucial point to 
observe here is that classical trajectories initiated in one of the regular regions cannot evolve into the other regular 
region. With increasing energy the classically disjoint regular regions move further apart and almost vanish near the 
dissociation energy. The result presented in Fig. [l] in fact closely mirrors the topological changes, shown in Fig. [2J in 
the phase space. Thus, in Fig. [l]a sequence of eigenstates with very small splittings begins right around the energy 
at which the symmetric periodic orbit becomes unstable. The unsymmetric stretch, however, becomes unstable at a 
higher energy E sa 8.3 and one can observe in Fig. [T] another sequence that seemingly begins near this point. Note 
that at higher energies it is not easy to identify any more sequences, but small splittings are still observed. The 
important thing to note here is that not only do the splittings but the individual eigenstates also correlate with the 
changes in the phase space. 

In Fig. [3]we show a set of four eigenstates to illustrate an important point - as soon as the 1:1 resonance manifests 
itself in the phase space, the tunneling doublets start to form and an integrable normal form approximation is 
insufficient to account for the splittings. Note that the phase space is mostly regular at the energies of interest. A 
much more detailed analysis can be found in the paper by Farrelly and Uzer[49^. We begin by noting [50] that the pair 
of eigenstates (counting from the zero-point) 35 and 36 are split by about AE sa 6.41 x 10~^ whereas the pair 37 and 
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FIG. 2. Evolution of the classical phase space with increasing energy shown as {u,pu) Poincare surface of sections. Note that 
large scale chaos appears for _E > 11. The formation of a separatrix and two classically disjoint regular regions can be seen at 
E = 9 due to the symmetric mode periodic orbit becoming unstable. The regular regions almost vanish near the dissociation 
energy. In the bottom row some of the chaotic orbits have been suppressed, for clarity of the figure, by showing them in gray. 



38 are separated by about 2.09 x 10"^. In both cases the splitting is smaller than the fundamental frequencies, which 
are of order unity. Note that the latter pair of states appears to be a part of the sequence in Fig. [ijthat starts right 
after the bifurcation of the symmetric stretch periodic orbit. Since the original mode frequencies are nonresonant, 
one can obtain a normal form approximation [15] to the original Hamiltonian and see if the obtained splittings can be 
explained satisfactorily. In other words, the observed AE are coming from the perturbation that couples both the 
modes and in such a case there is no reason for classifying them as tunneling doublets. However, from the coordinate 
space representations of the states shown in Fig. [3] one observes that there is an important difference between the 
two pairs of states. The pair (35, 36) seems to have a perturbed nodal structure and hence one can approximately 
assign the states using the zeroth-order quantum numbers (ns,nu). Inspecting the figure leads to the assignment 
(1, 6) and (0, 7) for states 35 and 36 respectively. If the above arguments are correct then the splitting AE'ss^sg should 
be obtainable from the normal form Hamiltonian. Since the theory of normal forms is explained in detail in several 
textbooks'51", I will provide a brief derivation below. Begin by using the unperturbed harmonic action-angle variables 
(I,</)) defined via the canonical transformation (similar set for s mode): 
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to express the original Hamiltonian as 
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where A = \^/2/{u]u\/^)- Since Hi is purely oscillatory the angle average 
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and thus to 0(A) the normal form Hamiltonian can be identified with the zeroth-order Hq(I) above. The first nontrivial 
correction arises at O(A^) and can be obtained using the 0(A) generating function 



Wi(I,c 
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where iJi^k(I) are the coefficients in the Fourier expansion of the oscillatory part of _ffi(I, 0). One now obtains the 
O(A^) correction to the Hamiltonian as 



K2 = -{Wi,Hi} 



(7) 



with the bar denoting angle averaging of the Poisson bracket involving Wi and Hi . Performing the calculations the 
normal form at O(A^) is obtained as 
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The primitive Bohr-Sommerfeld quantization Ij — > [rij + l/2)Ti yields the quantum eigenvalues perturbatively to 
0{\^). The procedure can be repeated to obtain the normal form Hamiltonian at higher orders. For instance, 
Farrelly and Uzer have[49] computed the normal form out to 0{X^^) and used Pade resummation techniques to 
improve in cases when the zeroth-order frequencies are near-resonant. For our qualitative discussions, the 0{\^) 
normal form is sufficient. 

Notice that the normal form Hamiltonian above is integrable since it is ignorable in the angle variables. Indeed, using 
the normal form and the approximate assignments of the eigenstates shown in Fig. ^ one finds A_E35^36 ~ 7.4 x 10"^ 
which is in fair agreement with the actual numerical value. Thus, a large part of the splitting Ai^as^ag can be 
explained classically. On the other hand, although the states 37 and 38 seem to have an identifiable nodal structure 
(cf. Fig. [3]), it is clear that they are significantly perturbed. Nevertheless, persisting with an approach based on 
counting the nodes, states 37 and 38 can be assigned as (8, 0) and (7, 1) respectively. Using the normal form one 
estimates Aii^37^38 « 5.12 x 10~^ and this is about a factor of two larger than the numerically computed value. Thus, 
the splitting in this case is not accounted for solely by classical considerations. One might argue that a higher order 
normal form might lead to better agreement, but the phase space Husimi distributions of the eigenstates shown in 
Fig. I3] suggests otherwise. The Husimi disributions, when compared to the classical phase spaces shown in Fig. [21 
clearly show that the pair (37,38) are localized in the newly created 1:1 resonance zone. Hence, this pair of states 
is directly influenced by the nonlinear resonance and the splitting between them cannot be accurately described by 
the normal form Hamiltonian. Indeed, as discussed by Farrelly and Uzer [33], in this instance one needs to consider 
a resonant Hamiltonian which explicitly takes the 1:1 resonance into account. This provides a clear link between 
dynamical tunneling, appearance of closely spaced doublets and creation of new, in this case a nonlinear resonance, 
phase space structures. The choice of states in Fig. [3] is different from what is usually shown as the standard example 
for dynamical tunneling pairs. However, in the discussion above the states were chosen intentionally with the purpose 
of illustrating the onset of near-degeneracy due to the formation of a nonlinear resonance. In a typical situation 
involving near-integrable phase spaces there are several such resonances ranging from low orders to fairly high orders. 
The importance of a specific resonance depends sensitively on the effective value of the Planck's constant [521 I53j. 
Indeed, detailed studies 54-56 have shown that excellent agreement with numerically computed splittings can be 
obtained if proper care is taken to include the various resonances. Clear and striking examples in this context can be 
found in the contributions by Schlagheck et al. and Backer et al. in this volume. 

To finish the discussion of the Davis-Heller system I show an example which involves states forming tunneling 
pairs that are fairly complicated both in terms of their coordinate space representations as well as their phase space 
Husimi distributions. The example involves three states 102, 103, and 104 around E w 14.4, which is rather close 
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FIG. 3. Four eigenstates close to an energy where the 1:1 resonance is just starting to appear. The upper panel shows the 
(m, s) coordinate space representations and the lower panels show the corresponding Husimi distributions in the {u,Pu) surface 
of section of the phase space. The first two states {E ~ 8.4) can be assigned approximate zeroth-order quantum numbers 
{us, Uu) whereas the last two states {E ~ 9.0) show perturbed nodal features. Clear difference in the phase space nature of the 
eigenstates can be seen. See text for details. 



to the dissociation energy. In the original workTT Davis and HcUcr noted that the splittings seem to increase by 
an order of magnitude in this high energy region. Notably, they commented that "perhaps the degree of irregularity 
between the regular regions plays some parf . In Fig. HI the variation of energy levels with the coupling parameter 
A is shown. One can immediately see that the three states are right at the center of an avoided crossing. The 
coordinate space representations of the states shows extensive mixing for states 102 and 104 while the state 103 
seems to be cleaner. The Husimi distributions for the respective states conveys the same message. Comparing to 
the phase space sections shown in Fig. [2] it is clear that the Husimis for state 102 and 104 seem to be ignoring the 
classical regular-chaotic divison - a clear indication of the quantum nature of the mixing. Although it is not possible 
to strictly assign these states as chaotic, a closer inspection does show substantial Husimi contribution in the border 
between the regular and chaotic regions. Interestingly, the pairwise splitting between the states is nearly the same 
but localized linear combinations of any two states exhibits two-level dynamics. Thus, the situation here is not of the 
generic [15] chaos-assisted tunneling one wherein one of the states is chaotic and interacts with the other two regular 
states. Nevertheless, linear combinations of the three states reveal (not shown here) that they are mixed with each 
other. A look at the coordinate and phase space representations of states 102 and 104 in Fig. [3] reveals that a different 
kind of state is causing the three-way interaction. It turns out that in the Davis-Heller model there is another class 
of symmetry-related pairs that appear when the unsymmetric mode becomes unstable. In the original work they 
were refered to as "circulating" states which displayed much larger splitting then the so-called local mode pairs. The 
case shown in Fig. [4] involves one of the circulating pairs interacting with the usual local mode doublet leading to 
the complicated three-way interaction. The subtle nature of this interaction is evident from the coordinate space 
representation of state 103 which exhibits a broken symmetry. Interestingly, and as far as I can tell, there has been 
very little understanding of such three-state interactions in the Davis- Heller system and further studies are needed to 
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FIG. 4. An avoided crossing in the high energy region {E ~ 14.4) of the Davis-Heller system involving three states. As in the 
previous figure, the {u, s) and phase space Husimi representations are also shown. The first and the last states appear to be 
strongly mixed states. The middle state seems relatively cleaner. Does chaos play a role in this case? Detailed discussions in 
the text. 



shed some light on the phase space nature of the relevant eigenstates. 



III. DYNAMICAL TUNNELING AND CONTROL: TWO EXAMPLES 



In the previous section the intimate connection between dynamical tunneling and phase space structures was 
introduced. The importance of a nonlinear resonance and a hint of the role played by the chaos (cf. Fig. [4]) is 
evident. Several other contributions in this volume discuss the importance of various phase space strucures using 
different models, both continuous Hamiltonians and discrete maps. In the molecular context, various mode-mode 
resonances play a critical role in the process of intramolecular vibrational energy redistribution (IVR)[57j. The 
phenomenon of IVR is at the heart of chemical reaction dynamics and it is now well established that molecules at 
high levels of excitation display all the richness, complexity and subtelty that is expected from nonlinear dynamics 
of multidimensional systems[58j. In this context, dynamical tunneling is an important agent [43l [59] of IVR and state 
mixing for a certain class of initial states (akin to the so called NOON states[SD]) which are typically prepared by the 
experiments. The importance of the anharmonic resonances to IVR and the regimes wherein dynamical tunneling, 
mediated by these resonances, is expected to be crucial is described in some detail in this volume by Leitner. The 
phase space perspective on Leitner's viewpoint can be found in a recent review[43] (see also Heller's contribution in the 
present volume) . In this regard it is interesting to note that there has been a renaissance of sorts in chemical dynamics 
with researchers critically examining the validity of the two pillars of reaction rate theory[61J - transition state theory 
(TST) and the Rice-Ramsperger-Kassel-Marcus (RRKM) theory. Since both theories have classical dynamics at their 
foundation, advances in our understanding of nonlinear dynamics and continuing efforts to characterize the phase 
space structure of systems with three or more degrees of freedom are beginning to yield crucial mechanistic insights 
into the dvnamics |62l I63j . At the same time, rapid advances in experimental techniques and theoretical understanding 
of the reaction mechanisms has led researchers to focus on the issue of controlling the dynamics of molecules. What 
implications might dynamical tunneling have on our efforts to control the atomic and molecular dynamics? In the rest 
of this article I focus on this issue and use two seemingly simple and well studied systems as examples to highlight the 
role of dynamical tunneling in the context of coherent control. Both examples are in the context of periodically driven 



systems and I refer the reader to the work of Flatte and Holthaus[Bl] for an exposition of the close quantum-classical 
correspondence in such systems. 

A. Driven quartic double well: chaos-assisted tunneling 

Historically, an early indication that dynamical tunneling could be sensitive to the chaos in the underlying phase 
space came from the study of strongly driven double well potential by Lin and Ballentine|17j. The model Hamiltonian 
in this case can be written down as 

H{x,p; t) — Hq{x,p) + Xix cos{uj pt) , (10) 

with ujp being the frequency of the monochromatic field (henceforth refered to as the driving field) and the unperturbed 
part 

Ho{x,p) = ^p^+Bx^-Dx^ (11) 

is the Hamiltonian corresponding to a double well potential with two symmetric minima at a: = ±(Z3/2i?)^'^ and a 
maximum at x = 0. Following the original work[17j. the parameters of the unperturbed system (assuming atomic 
units) are taken tobeM = l,i? = 0.5, and D — 10 for which the potential has a barrier height Vb = 50 and supports 
about eight tunneling doublets. As is well known, in the absence of the driving field, a wavepacket prepared in the 
left well can coherently tunnel into the right well with the time scale for tunneling being inversely proportional to 
the tunnel splitting. This unperturbed scenario, however, is significantly altered in the presence of a strong driving 
field. In the presence of a strong field the phase space of the system exhibits large scale chaos coexisting with two 
symmetry-related regular regions. Lin and Ballentine observed that a coherent state localized in one of the regular 
region tunnels to the other symmetry-related regular region on timescales which are orders of magnitude smaller than 
in the unperturbed case. It was suspected that the extensive chaos in the system might be assisting the tunneling 
process. 

In order to illustrate the tunneling process Fig. [5^ shows the stroboscopic surface of section for the case of strong 
driving with Ai = 10. Note the extensive chaos and the two regular islands (left and right) in the phase space. A 
coherent state \z) = \xo,po) is placed in the center of the left island and time evolved using the Floquet approach 
which is ideally suited for time-periodic driven systems. In this instance one is interested in the time at which the 
coherent state localized on the left tunnels over to the regular region on the right. In order to obtain this information 
it is necessary to compute the survival probability of the initial coherent state. Briefly, Floquet states {|Xn)} are 
eigenstates of the Hermitian operator H — ifid/dt and form a complete orthonormal basis. An arbitrary time-evolved 
state \z{t)) can be expressed as 

|z(t))=^A„e-^^"*|x„(i)), (12) 

n 

with En being the quasienergy associated with the Floquet state \xn)- The expansion coefficients An are independent 
of time and given by 

An - (x„(0)|z(0)), (13) 

yielding the expansion 

\z{t)) =Y.e-^^-'\Xn{t)){Xnmzm- (14) 

n 

Measuring time in units of field period (T/) and owing to the periodicity of the Floquet states, \Xn{'t)) = \Xn{'t + Tf))^ 
the above equation simplifies to 

\z{t)) = Y,e-^^-^\xnm{xnmzm ^u{T)\zm, m 



with T — kTf and integer k. 
The time evolution operator 



Uir) = J2e-^^-^\Xnm{Xnm, (16) 



10 



20 


1 1 1 1 1 , 1 1 1 1 1 


TO 


//^''^X " 


P 




-to 


- 


-ao 


: 1 . 1 , 1 , 1 . 1 . - 




10 



10' 




FIG. 5. (a) Phase space for the driven double weU system with field strength Ai — 10. A coherent state localized in the 
left regular island tunnels (indicated by a green arrow) over to the right regular island. The chaotic regions (gray) have been 
suppressed for clarity, (b) Monitoring the survival probability of the initial state to determine the timescale of tunneling 
(~ lOOr in this case). Snapshots of the evolving Husimi distributions are also shown at specific intervals, with green (gray) 
indicating maxima of the distributions, (c) The decay time determined as in (b) for a variety of field strengths with the initial 
state localized in the left island. Note the fluctuations in the decay time over several orders of magnitude despite very similar 
nature of the phase spaces over the entire range of the driving field strength. The arrows highlight some of the plateau regions 
which are crucial for bichromatic control. 



is determined by successive application of the one-period time evolution operator U{Tf, 0) i.e., 

U{kTf,0) = [U{TfM''- (17) 

This allows us to express the survival probability of the initial coherent state in terms of Floquet states as 

5(r)^|(z(0)|z(r))p 



E< 



-iE„T 



^PznPz 



:^(0)|Xn(0))(Xn(0)|z(0)) 



i(E„-_E„)r 



(18) 



where the overlap intensities are denoted by Pzn = l('2^(0)|Xn(0))p. In order to determine the "lifetime" (T) of the 
initial coherent state, we monitor the time at which S{t) goes to zero (minimum) at the first instance. In other words. 
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FIG. 6. Decay time plot for ujf = ijJp = 6.07 (black) and for u)p + 0.05 (blue) and iOp — 0.05 (red) as a function of Ai. 
The local phase structure in the vicinity of the initial coherent state and the survival probability for fixed Ai = 6.4 (indicated 
by green line in (a)) corresponding to red shifted, fundamental and blue shifted field frequency is shown in panels (b) and (c) 
respectively. The 6:1 field-matter resonance in case of ujp — 0.05 is clearly visible and correlates with the very short decay 
time as compared to the other two cases. Note the different time axis scale in the first survival probability plot. 

it is the time when the coherent state leaves its initial position i. e, the left regular island of the phase space for the 
first time. In Fig. [sb the S{t) for the initial state of interest is shown along with the snapshots of Husimi distribution 
at specific times. Clearly, the initial state tunnels in about lOOT/, which suggests that chaos assists the dynamical 
tunneling process. However, the issue is subtle and highlighted in Fig. [5J: which shows the decay time plot for a range 
of driving field strengths for an intial state localized in the left regular island. It is important to note that the gross 
features of the phase space are quite similar over the entire range. However, Fig. [Sb exhibits strong fluctuations over 
several orders of magnitude and this implies that a direct association of the decay time with the extent of chaos in 
the phase space is not entirely correct. 



The above discussion and results summarized in Fig. [5] bring up the following key question. What is the mechanism 
by which the initial state \z) decays out of the regular region? In turn, this is precisely the question that modern 
theories of dynamical tunneling strive to answer. According to the theory of resonance-assisted tunneling [Ml [251 I52j . 
the mechanism is possibily one wherein | z) couples to the chaotic sea via one or several nonlinear resonances provided 
certain conditions are satisfied. Specifically, the local structure of the phase space surrounding the regular region is 
expected to play a critical role. The theoretical underpinnings of RAT along with several illuminating examples can 
be found in the contribution by Schlagheck et al. and here we suggest a simple numerical example which points to 
the importance of the local phase space structure around \z). 

Preliminary evidence for the role of field-matter nonlinear resonances in controlling the decay of \z) is given in 
Fig. [6j which shows the effect of changing the driving field frequency cop on the local phase space structures and the 
decay times. From Fig. [6] it is apparent that detuning ujf by ±0.05 leads to a significant change in the local phase 
space structure and the decay time of the initially localized state. In Fig. pfb), corresponding to the field frequency 
Ldp — 6.02, a prominent 6:1 field-matter resonance is observed. It is plausible that the the decay time is only a few 
hundred field periods in this case due to assistance from the nonlinear resonance. However, the decay time increases 
for top = 6.07 and becomes even larger by an order of magnitude for ujp = 6.12, due to absence of the 6:1 resonance. 
This indicates that decay dynamics is highly sensitive to the changes in local phase space structure of the left regular 
island. Therefore, taking into account the relatively large order resonances, since ?i = 1, is unavoidable in order to 
understand the decay time plot in Fig. [5] and for smaller values of the effective Planck constant one expects a more 
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FIG. 7. Stroboscopic surface of sections for the bichromatically driven double well system with the control field strength fixed 
at A2 = 0.7 and primary field strength (a) Ai — 9.2 and (b) Ai = 10.0. As in the earlier plots, the chaotic regions have been 
suppressed (in gray) for clarity. Note the breaking of the symmetry in both cases. However, significant suppression of the deacy 
of the state localized in the left regular island happens in case (b) only. 

complicated behavior. Note, however that the very high order island chain visible in the last case in Fig.[6]is unable to 
assist the decay. This is where we believe that an extensive 7i-scaling study will help in gaining a deeper understanding 
of the decay mechanism. Such an extensive calculation can be found, for example, in the recent work [54] by Mouchet, 
Eltschka, and Schlagheck on the driven pendulum. There is sufRcient evidence [54] in the driven pendulum system 
for a mechanism in which nonlinear resonances play a central role in coupling initial states localized in regular phase 
space regions to the chaotic sea. One might be able to provide a clear qualitative and quantitative explanation for 
the results in Fig. [5t based on the recent advances. 



CAT spoils bichromatic control 



This brings us to the second important issue - what is the role of chaos in this dynamical tunneling process? An 
earlier work[HS] by Utermann, Dittrich and Hanggi on the driven double well system showed that there is indeed a 
strong correlation between the splittings of the Floquet states and the overlaps of their Husimi distributions with the 
chaotic regions in the phase space. However, a different perspective yields clear insights into the role of chaos with 
potential implications for coherent control. In order to highlight this perspective I start with a simple question: is it 
possible to control the decay of the localized initial state with an appropriate choice of a control field? It is important to 
note that the Lin-Ballentine system parameters imply that one is dealing with a multilevel control scenario and there 
has been a lot of activity over the last few years to formulate control schemes involving multiple levels in both atomic 
and molecular systems. The driven double well system has been a particular favorite in this regard, more so in recent 
times due to increased focus on the physics of trapped Bose-Einstein condensates [HS]. More specifically, several studies 
have explored the possibility of controlling various atomic and molecular phenomenon using bichromatic fields with 
the relative phase between the fields providing an additional control parameter. In the present context, for example, 
Sangouard et al. exploited the physics of adiabatic passage to show that an appropriate combination of (a;F,2cjj?) 
field leads to suppression of tunneling in the driven double well model(67|. The choice of relative phase between the 
two fields allowed them to localize the initial state in one or the other well. However, the parameter regimes in the 
work by Sangouard et al. correspond to the underlying phase space being near-integrable and hence a minimal role 
of the chaotic sea. 
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More relevant to the mixed regular-chaotic phase space case presented in Fig.lslis an earlier work'SS' by Farrelly and 
Milligan wherein it was demonstrated that one can suppress the tunneling dynamics in a driven double well system 



using a {ujf, 2w^) bichromatic field. In other words the original Hamiltonian of Eq. 10 is modified as following 



H{x,p;t) — H(){x,p) + Xixcos{ujFt) + \2XCOs{2ujpt + cp), (19) 

with Ho{x,p) being the same as in Eq. [TT|and the additional 2wi?-field is taken to be the control field. Moreover, 
modulating the turn-on time of the control field can trap the wavepacket in the left or right well of the double well 
potential. Hence, it was argued J68] that the tunneling dynamics in a driven double well can be controlled at will for 
specific choices of the control field parameters (A2, (f>)- Note that in the presence of control field i.e., A2 y^ with the 
relative phase = the Hamiltonian in Eq. [19] transforms under symmetry operations as 

H I —X, —p; t -\ I = Ho{x,p) — x[Xi cos{ujpt + tt) + A2 cos{2uJFt + 27r)] 

= Hq{x,p) — x[— Ai coa{ujFt) + A2 cos(2wi?t)] (20) 

^H{x,p;t). 

Similarly, except at = 7''/2, the discrete symmetry of the Hamiltonian is broken under the influence of the additional 
2a;i;'-field. Farrelly and Milligan thus argued[S5] that the control field with strength smaller than the driving field will 
lead to localization due to the breaking of the generalized symmetry of the Hamiltonian and the Floquet states. The 
impact of a small symmetry breaking control field with (A2, 0) = (0.7, 0) can be clearly seen in Fig. [t] in terms of the 
changes in the classical phase space structures. However, there is a subtlety which is not obvious upon inspecting 
the phase spaces shown in Fig. [7] for two different but close values of the driving field strength. In case of Ai = 9.2 
corresponding to Fig. l7[a), the control field is unable to suppress the decay of the initial state \z) localized in the left 
regular region. On the other hand. Fig. WCb) corresponds to Ai = 10 and computations show that the control field is 
able to suppress the decay of \z) to an appreciable extent. Thus, although in both cases the control field breaks the 
symmetry and the resulting phase spaces show very similar structures, the extent of control exerted by the 2wf^-field 
is drastically different. 

The discussions above and the results summarized in Fig. [t] and Fig. [6] clearly indicate that the decay of \z) out of 
the left regular region can be very different and deserves to be understood in greater detail. In fact, computations 
show that such cases of complete lack of control are present for other values of Ai as well. This is confirmed by 
inspecting Fig. [8] which shows the control landscape for the bichromatically driven double well in the specific case of 
(j) — 0. Other choices for (p also show similarly convoluted landscapes. There are several ways of presenting a control 
landscape and in Fig. [S] the time-smoothed survival probability (cf. Eq. 18) 



(S)^- f dr'Sir') (21) 

associated with the initial state \z) is used to map the landscape as a function of the field strengths (Ai, A2). Note 
that the choice of {S) to represent the landscape is made for convenience; the decay time is a better choice which 
requires considerable effort but the gross qualitative features of the control landscape do not change upon using {S). 
Large (small) values of {S) indicate that the decay dynamics is suppressed (enhanced). It is clear from Fig. [S] that the 
landscape comprises of regions of control interspersed with regions exhibiting lack of control. Such highly convoluted 
features of control landscape are a consequence of the simple bichromatic choice for the control field and the nonlinear 
nature of the corresponding classical dynamics. From a control point of view, there are regions on the landscape for 
which a monotonic increase of A2 leads to increasing control. Interestingly, the lone example of control illustrated in 
Farrelly and Milligan's paper 68^ happens to be located on one of the prominent hills on the control landscape (shown 
as a green dot in Fig. [s]). A striking feature that can be seen in Fig. [s] is the deep valley around Ai — 9.2 which 
signals an almost complete lack of control even for significantly large strengths of the 2wi?-field. The valleys in Fig. [8] 
correspond precisely to the plateaus seen in the decay time plot shown in Fig. p^ for A2 = (red arrows). It is crucial 
to note that this "wall of no control" is robust even upon varying the relative phase (j) between the driving and the 
control field. 

Insights into the lack of control for driving field strength Ai = 9.2 (and other values as well which are not discussed 
here) can be obtained by studying the variation of the Floquet quasienergies with A2 , the control field strength. In 
Fig. [9] the results of such a computation are shown and it is immediately clear that even in the presence of 2wi?-control 
field six states contribute to the decay of \z) over the entire range of A2. This is further confirmed in Fig. [91(b) where 
the plot of the overlap intensities shows multiple Floquet states participating nearly equally in the decay dynamics 
of the initial state. The final clue comes from inspecting the Husimi distributions shown in Fig. [9| highlighting the 
phase space delocalized nature of some of the participating states. Despite the symmetry of the tunneling doublets 
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FIG. 8. Time-smoothed survival probability representation of the control landscape for (f} — as a, function of the field 
parameters (Ai,A2). The initial state in every case is localized in the left regular region in the classical phase space. Notice 
the convoluted form of the landscape with the regions of low probability indicating little to no control (red thick arrows). The 
green lines (thin arrows) indicate regions where a high degree of control can be achieved. 



being broken due to the bichromatic field, two or more of the participating states are extensively delocalized in the 
chaotic reigons of the phase space. Moreover, the participation of the chaotic Floquet states persists even for larger 
values of A2. Hence, using the symmetry breaking property of the 2aji?-field for control purposes is not very effective 
when chaotic states are participating in the dynamics. Therefore, the lack of control, signaled by plateaus in Fig. [^ 
and the valleys in Fig. [^ is due to the dominant participation by chaotic states i.e., chaos-assisted tunneling. The 
plateaus arise due to the fact that the coupling between the localized states and the delocalized states vary very little 
with increasing control field strength A2 - something that is evident from the Floquet level motions shown in Fig. [9] 
and established earlier by Tomsovic and UUmo in their seminal work on chaos-assisted tunneling in coupled quartic 
oscillators HH]. It is important to note that for A2 = the chaotic states, as opposed to the regular states, do not have 
a definite parity. Consequently, the presence of the 2wi;--field does not have a major influence on the chaotic states. 
Thus, if one or more chaotic states are already influencing the dynamics of \z) at A2 = then the bichromatic control 
is expected to be difficult. An earlier study^69, by Latka et al. on the bichromatically driven pendulum system also 
suggested that the ability to control the dynamics is strongly linked to the existence of chaotic states. 

The model problem in this section and the results point to a direct role of chaos-assisted pathways in the failure 
of an attempt to bichromatically control the dynamics. However, it is not yet clear if control strategies involving 
more general fields would exhibit similar characteristics. There is some evidence in the literature which indicates that 
quantum optimal control landscapes might be highly convoluted if the underlying classical phase space exhibits large 
scale chaos[70]. Nevertheless, further studies need to be done and the resulting insights are expected to be crucial in 
any effort to control the dynamics of multilevel systems. 



B. Driven Morse oscillator: resonance-assisted tunneling 



The driven Morse oscillator system has served as a paradigm model for understanding the dissociation dynamics 
of diatomic molecules. Studies spanning nearly three decades have explored the physics of this system in exquisite 
detail. Consequently, a great deal is known about the mechanism of dissociation both from the quantum and classical 
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FIG. 9. (a) Variation of the Floquet quasienergies with A2. The primary driving field strength is fixed at Ai — 9.2. Six states 
that have appreciable overlap with \z) are highlighted by circles, (b) Overlap intensity p^n for Ai — 9.2 and A2 = 2.1 indicates 
multilevel interactions involving the states shown in (a). Husimi distribution function of the Floquet states regulating the decay 
of \z) are also shown. Notice that the nature of delocalized states (magenta) does not change much with A2. 



perspectives. Indeed, the focus of researchers nowadays is to control, either suppress or enhance, the dissociation 
dynamics and various suggestions have been put forward. In addition, one hopes that the ability to control a single 
vibrational mode dynamics can lead to a better understanding of the complications that arise in the case of polyatomic 
molecular systems wherein several vibrational modes are coupled at the energies of interest. 

Several important insights have originated from classical-quantum correspondence studies which have established 
that molecular dissociation, in analogy to multiphoton ionization of atoms, occurs due to the system gaining energy by 
diffusing through the chaotic regions of the phase space. For example, an important experimental study by Dietrich and 
Corkum has shown 71j, amongst other things, the validity of the chaotic dissociation mechanism. Thus, the formation 
of the chaotic regions due to the overlap [75] of nonlinear resonances (field-matter), hierarchical structures [^ near 
the regular-chaotic borders acting as partial barriers, and their effects on quantum transport I74j have been studied 
in a series of elegant papers [75U77J . In the context of this present article an interesting question is as follows. Since a 
detailed mechanistic understanding of the role of various phase space structures of the driven Morse system is known, 
is it possible to design local phase space barriers to effect control over the dissociation dynamics? In particular, the 
central question here is whether the local phase space barriers are also able to suppress the quantum dissociation 
dynamics. There is an obvious connection between the above question and the theme of this volume - quantum 
mechanics can "shortcircuit" the classical phase space barriers due to the phenomenon of dynamical tunneling. Thus, 
such phase space barriers might be very effective in controlling the classical dissociation dynamics but might fail 
completely when it comes to controlling the quantum dissociation dynamics. There is a catch here, however, since 
there is also the possibility that the cantori barrier in the classical phase space may be even more restrictive in the 
quantum case. Thus, creation or existence of a phase space barrier invariably leads to subtle competition between 
classical transport and quantum dynamical tunneling through the barrier. As expected, the delicate balance between 
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FIG. 10. A plot of Morse potential along with the bound states supported by the potential well. A schematic view of energy 
level diagram for Morse eigenstates under study in this chapter is shown as inset. The laser field of intensity Ai — 0.0287 and 
frequency up = 0.0178 connects state n = 10 and n — 14 via a two-photon resonance transition. The location of cantori in 
presence of the field with top = 0.0178 are indicated as dotted curves in the figure. 

the classical and quantum mechanism is determined by the effective Planck constant of the system of interest. An 
earlier detailed review[32] by Radons, Geisel and Rubner is highly reccomended for a nice introduction to the subject of 
classical-quantum correspondence perspective on phase space transport through Kolmogoroff-Arnol'd-Moser (KAM) 
and cantori barriers. In the driven Morse oscillator case, Brown and Wyatt showed [75] that the cantori barriers 
do leave their imprint on the quantum dissociation dynamics and act as even stronger barriers as compared to the 
classical system. Maitra and Heller in their studv[33l on transport through cantori in the whisker map have clearly 
highlighted the classical versus quantum competition. 



From the above discussion it is apparent that any approach to control the dissociation dynamics by recreating local 
phase space barriers will face the subtle classical-quantum competition. In fact, it is tempting to think that every 
quantum control algorithm works by creating local phase space dynamical barriers and the efficiency of the control is 
decided by the classical-quantum competition. However, at this point of time there is very little work towards making 
such a connection and the above statement is, at best, a conjecture. For the purpose of this article, I turn to the 
driven Morse oscillator system to provide an example for the importance of resonance- assisted tunneling in controlling 
the dissociation dynamics. 

The model system is inspired from the early work [75] by Wyatt and Brown (see also the work by Breuer and 
Holthaus'78') and the Hamiltonian can be written as 



H{x,p; t) = Hq(x,p) — /i(x)Ai cos{ujpt), 
in the dipole approximation. The zeroth-order Hamiltonian 



Ho{x,p) 
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2M 



P 



Do[l 



-a(2;— a;e)i 



(22) 



(23) 
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represents the Morse oscillator modeling the anharmonic vibrations of a diatomic molecule. In the above, ii{x) is the 
dipole moment function, Ai is the strength of the laser field, up is the driving field frequency, and 

M= -^-, (24) 

mi + 1712 



is the reduced mass of the diatomic molecule with rrii and m2 being the two atomic masses. In Eq. 23 Dq is the 
dissociation energy, a is the range of the potential and Xe is the equilibrium bond length of the molecule. 

Rather than attempting to provide a general account as to how RAT might interfere with the process of control, 
I feel that it is best to illustrate with a realistic molecular example. Hopefully, the generality of the arguments 
will become apparent later on. For the present purpose I choose the diatomic molecule hydrogen fluoride (HF) as 
the specific example. Any diatomic molecule could have been chosen but HF is studied here due to the fact that 
Brown and Wyatt have already discussed the role of cantori barriers to the dissociation dynamics in some detail. The 
Morse oscillator parameters for hydrogen fluoride are Dq — 0.225, a — 1.174, Xe = 1.7329 and M = 1744.59. These 
parameters correspond to ground electronic state of the HF molecule supporting Nb = 24 bound states. Note that 
atomic units are used for both the molecular and field parameters with time being measured in units of the field 
period Tp = 2tt/ujp. The only difference between the present work and that of Brown and Wyatt has to do with the 
field-matter coupling. Brown and Wyatt use the dipole function 

fiix) ^ Axe-^'^', (25) 

with A = 0.4541 and B = 0.0064, obtained from ab-initio data on HF. Here a linear approximation for j-i{x) 

H{x) fa fi{Xe)+ i-K-j {x~Xe) 

= ^i{Xe) + di{x - Xe), (26) 

is employed with di « 0.33 in case of HF. There are quantitative differences in the dissociation probabilities due to the 
linearization approximation but the main qualitative features remain intact despite the linearization approximation. 



In Fig. 10 the Morse potential for HF is shown along with a summary of the key features such as the energy region 
of interest, the quantum state whose dissociation is to be controlled, and the classical phase space structures that 
might play an important role in the dissociation dynamics. The driving field parameters are chosen as {lup,Xi) — 
(0.0178,0.0287), same as in the earlier workTF, and the field strength corresponds to about 30 TW/cm^ = 30 x lO^'^ 
W/cni'^. As shown in Fig. 10 (inset), the focus is on understanding and controlling the dissociation dynamics of the 



n = 10 excited Morse oscillator eigenstate. There are several reasons for such a choice and I mention two of the 
most important reasons. Firstly, the earlier studyfTS] has established that n = 10 of HF happens to be in an energy 
regime wherein two specific cantori barriers in the classical phase space affect the dissociation dynamics. Moreover, 
for driving frequency ujp = 0.0178, two of the zeroth-order eigenstates n = 10 and n — 14 have unperturbed energies 
such that Ei4 — Eiq w 2hujp and hence corresponds to a two-photon resonant situation. Secondly, for a field strength 
of 30 TW/cm^, the dissociation probability of the ground vibrational state is negligible. Far stronger field strengths 
are required to dissociate the n = state and ionization process starts to compete with the dissociation at such high 
intensities. Thus, in order to to illustrate the role of phase space barriers in the dissociation dynamics without such 
additional complications, the specific initial state n = 10 is chosen. Incidentally, such a scenario is quite feasible since 
a suitably chirped laser field can populate the n = 10 state very efficiently from the initial ground state n — and 
one imagines coming in with the monochromatic laser to dissociate the molecule. 

1. Nature of the classical phase space 

The monochromatically driven Morse system studied here has a dimensionality such that one can conveniently visu- 
alize the phase space in the original cartesian {x,p) variables. However, since the focus is on suppressing dissociation 
by creating robust KAM tori in the phase space, action-angle variables (J, 6) which are canonically conjugate to {x,p) 
are convenient and a natural representation to work with. The action- angle variables (J, 9) of the unperturbed Morse 
oscillator, appropriate for the bound regions, are given by [75] 



^-\/^(l-^/^^^) (27a) 

(27b) 



9 = ~sgn{p) cos 



'^-^ ^a(x-x,) 1 



E V-B 
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In the above equations, E = Hq/Dq < 1 denotes the dimensionless bound state energy, and sgn{p) = 1 for p > 0, 
sgn{p) = —1 for p < 0. In terms of the action-angle variables it is possible to express the cartesian {x,p) as follows 



X — Xe -\ In 

a 



1 + ^yEo{J)cose 



P 



(l-^o(J)) 

1 + ^/Eq{j) cose 



(28a) 
(28b) 



where Eq{J) = Hq{J)/Dq. Substituting for {x,p) in terms of (J, 6*), the unperturbed Morse oscillator Hamiltonian in 
Eq. 23 is transformed into 



ffo(J)=-o(J-^J^ 



(29) 



where wq = ^J^o^Dq/M is the harmonic frequency at the minimum. The zeroth-order nonlinear frequency is easily 
obtained as 



o„(.).fi..„(i-|^. 



(30) 



and with increasing excitation i.e., increasing action (quanta) J, the non-linear frequency ^o{J) decreases monoton- 
ically and eventually vanishes, signaling the onset of unbound dynamics leading to dissociation. 



The driven system can now be expressed in terms of the variables (J, ^) as 

1 + y/Eo{J) con 6 



H(J,e;t) ^Ho{J)- -In 
a 



In addition, since a; is a periodic function of 9, one has the Fourier expansion 



cos{ujFt). 



x = 2 



K)(J) + 5Z^"('^)cos(r 



n=l 



As a consequence the driven Hamiltonian in Eq. [31] can be written as 

HiJ,e;t)^HoiJ)-eviJ,e;t), 
where the matter-field interaction term is denoted as 



viJ,e;t) 



VoiJ) + J2M J) cosine) 



cos(aji?t). 



The Fourier coefficients Vo(J) and Vn{J) are known analytically and given by 



W) = ^ln 



Dn 



Do^-DoEoiJ) 



2{Do-Eo{J)) 



VniJ) 
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(31) 



(32) 
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(34a) 



(34b) 



The stroboscopic surface of section in the (J, e) variables is shown in Fig. [TT] and is a typical mixed regular-chaotic 
phase space. A few important points are worth noting at this stage. Firstly, the initial state of interest is located 
close to a cantorus with ojf '■ ^oiJ) = 1 + 7J^ with 7g = (1 + \/5)/2 being the golden ratio. The importance of 
this cantorus to the resulting dissociation dynamics of the n = 10 state was the central focus of the work by Brown 
and Wvatt [75] . In particular, the extensive stickiness around this region can be clearly seen and hence one expects 
nontrivial influence on the classical dissociation dynamics as well. Secondly, a prominent ujp:flo{J) = 2:1 nonlinear 
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FIG. 11. The phase space, as a stroboscopic surface of section, for the driven Morse system with laser field of intensity 
Ai = 0.0287 and frequency ojf = 0.0178 corresponding to a on-resonance situation. The 2-photon resonance is clearly seen in 
the phase space and marked as 2:1 resonance in dark green. The initial Morse eigenstate n = 10 (thick dashed line) is situated 
rather close to a cantorus (indicated), ujf/^o{J) ~ 1 + 7^^ with 7g = (1 + \/5)/2 being the golden ratio. The n = 10 state 
is connected to the n = 14 Morse eigenstate via the 2:1 resonance. Note that the states n — 10, 12, and 14 are symmetrically 
located about the resonance with n = 12 being localized in the resonance. 



resonance is also observed in the phase space and represents the classical analog of the quantum 2-photon resonance 
condition. Interestingly, the area of this resonance is about h and, therefore, can support one quantum state. It turns 
out that the Husimi density of the Morse state n = 12 is localized inside the 2:1 resonance island. In addition, the 
states n = 10 and n — lA are nearly symmetrically located about the 2:1 resonance. A way to see this is to use 
secular pertubation theory on the driven Morse Hamiltonian. One can show that in the vicinity of the 2:1 resonance 
(cf. Fig. 11) an effective pendulum Hamiltonian 



Zm2:l 

is obtained with A J = J — J2:i, m2:i = IDqJuj^ and V2:i(J) = eVi(J)/2. The resonant action 
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(35) 
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for the parameters used in this work. Using action values J == 10.5 (quantum state n — 10) and J' = 14.5 (quantum 
state v! = 14) the energy difference is calculated as 



\Ej - Ej. 



1 



(J- j')(J+J'-2J2:l) 



2m2:l 

3.2 X lO"** < Ej 



(37) 



In other words, the states n = 10 and n = 14 are nearly symmetrical with respect to the state n = 12, which is 
localized in the 2:1 resonance. Therefore, the nonzero coupling V2:i will efficiently connect the states n — 10 and 



20 

n = 14. Moreover, for the given parameters, using the definition of Vr.siJ) in terms of Fourier coefficient Vi, the 
strength of the resonance is estimated as T^2:i(</2:i) ~ 0.01464, clearly a fairly strong resonance. Consequently, the 
situation in Fig. [TT] is a perfect example where RAT can play a crucial role in the dissociation dynamics. Indeed 
quantum computations (not shown here) show that there is a Rabi-type cycling of the probabilities between the three 
Morse states. I now turn to the issue of selectively controlling the dissociation dynamics of the initial Morse eigenstate 
n = 10 by creating local barriers in the phase space shown in Fig. [Ill bearing in mind the possibility of quantum 
dynamical tunneling interfering with the control process. 

2. Creating a local phase space barrier 

If the classical mechanism of chaotic diffusion leading to dissociation holds in the quantum domain as well then a 
simple way of controlling the dissociation is to create a local phase space barrier between the state of interest and 
the chaotic region. In a recent work(50], Huang, Chandre and Uzer provided the theory for recreating local phase 
space barriers for time-dependent systems and showed that such barriers indeed suppress the ionization of a driven 
atomic system. However, Huang et al. were only concerned with the classical ionization process. Thus, potential 
complications due to dynamical tunneling were not addressed in their study. The driven Morse system studied here 
presents an ideal system to understand the interplay of quantum and classical dissociation mechanisms. In what 
follows, I provide a brief introduction to the methodology with an explicit expression for the classical control field 
needed to recreate an invariant KAM barrier, preferably an invariant torus with sufficiently irrational frequency ^r- 

To start with, the nonautonomous Hamiltonian is mapped into an autonomous one by considering (i(mod27r),_E) 
as an additional angle-action pair. Denoting the action and angle variables by A = {J,E) and 9 = {9,t), the original 



driven system Hamiltonian (see Eq. 33 1 can be expressed as 

HiA,e)=HoiA)-eV{A,e), (38) 

with V{A,6) = v{J,9;t). Note that for a fixed driving field strength Ai and the value of di corresponding to a 
diatomic molecule, e = Xidi is also fixed. Moreover, for physically meaningful values of di for most diatoms and 
typical field strengths far below the ionization threshold one always has e ^ 1. In the absence of the driving field 
(e = 0), the zeroth-order Hamiltonian is integrable and the phase space is foliated with invariant tori labeled by the 
action A corresponding to the frequency a; = dH^/dA = (I^qjI^f)- However, in the presence of the driving field 
(e ^ 0) the field-matter interaction renders the system nonintegrable with a mixed regular-chaotic phase space. More 
specifically, for field strengths near or above a critical value Ec one generally observes a large scale destruction of the 
field-free invariant tori leading to significant chaos and hence the onset of dissociation. The critical value tc itself is 
clearly dependent on the specific molecule and the initial state of interest. The aim of the local control method is to 
rebuild a nonresonant torus Aq = (Jo,0), k • a; ^ with integer k, which has been destroyed due to the interaction 
with the field. Assuming that the destruction of Aq is responsible for the significant dissociation observed for some 
initial state of interest, the hope is that locally recreating the Aq will suppress the dissociation i.e., Ao acts as a 
local barrier to dissociation. Ideally, one would like to recreate the local barrier by using a second field (appropriately 
called as the control field) which is much weaker and distinct from the primary driving field. 

Following Huang[80 et al. such a control field f{6) can be analytically derived and has the form 

f{e) = -H{Ao^d0Th{d),e), (39) 

where h{6) = H{Aq, 9) = J^k by^e'^^'^ and F being a linear operator defined by 

r^^) - E ^^""- (40) 

The classical control Hamiltonian can now be written down as 

H,{A,e) = H{A,e) + f{e) 

= HoiA)-eV{A,9) + f{e). (41) 

In case of the driven Morse system the control field can be obtained analytically and to leading order is given by 

H,iJ,9;t) ^ Ho{J) - eviJ,9;t) + e^ga{e,t), (42) 

where 

9ai9,t) = ^{deTbf - 2Voi{derb)cos{ujFt) - {derb)aJ,9;t), (43) 
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and it can be shown that 



V'oi — 



vn. = (-1)"- (^) i"°"^-;!:; , (44) 

\2aDoJ {oja + Qr)" 

00 

C( J, 0; i) = ^ Vni{Jr)[{cos{n6 + ixjpt) + cos(n6' - LUpi)] . 

71 = 1 

I skip the somewhat tedious derivation of the resuh above and refer to the original hterature[80 as well as a recent 
thesis [454 for details. Note that in the above fir is the frequency of the invariant torus that is to be recreated 
corresponding to the unperturbed action 

or) 

Jr = {uJo-^r) 2", (45) 

and to 0(e) is located at J{6) — Jr ~ edoTb, assuming the validity of the perturbative treatment. 



The leading order control field in Eq. 43 is typically weaker than the driving field and has been shown 45 to be 
quite effective in off-resonant cases in suppressing the classical dissociation. However, in order to study the effect of 
the control field on the quantum dissociation probabilities, it is necessary to make some simplifications. One of the 
main reasons for employing the simplified control fields via the procedure given below has to do with the fact that 
the classical action- angle variables do not have a direct quantum counterpart |81j. Essentially, the dominant Fourier 



modes Fk-^^k^ of Eq. 43 are identified and one performs the mapping 

Fki,k2cos{ki0 + k2LUFt) -^ A2(fci,fc2)cos(fc2WFi), (46) 

yielding the simplified control Hamiltonian 

He = H{,1 e- 1) + A*(a;)A2(fci, fcs) cos(fc2WFt). (47) 

If more than one dominant Fourier modes are present then they will appear as additional terms in Equation |47[ Note 
that the above simplified form is equivalent to assuming that the control field is polychromatic in nature, which need 
not be true in general. Nevertheless, a qualitative understanding of the role of the various Fourier modes towards 
local phase space control is still obtained. More importantly, and as shown next, in the on-resonant case of interest 
here, the simplified control field already suggests the central role played by RAT. 

3. RAT spoils local phase space control 

In Fig. [12] a summary of the efforts to control the dissociation of the initial n = 10 state is shown. Specifically, 
Fig. 12 ^a) and (b) show the phase spaces where two different KAM barriers, cjp/ilr = 1 -f 7"^ and ujp/flr = VS 



respectively are recreated. The control Hamiltonian in both instances, in cartesian variables, has the following form 

He = Ho{x,p) — \i^{x) cos{LLipt) + \21jl{x) cos{2ujpt), (48) 

with A2 ~ 0.01 = 3 TW/cm^^. Thus, as desired, the control field strengths are an order of magnitude smaller than 
the driving field strength. The simplified control Hamiltonian reflects the dominance of the F3__2 Fourier mode of the 
leading order control field in Eq. [43l and obtained using the method outlined before. Also, note that the control field 
comes with a relative phase (j) = tt with respect to the driving field. As expected, from the line of thinking presented 
in the previous section. Fig. ll2Fc) and (d) show that the KAM barriers indeed suppress the classical dissociation 
significantly. However, the quantum dissociation in both cases increases slightly! Clearly, the recreated KAM barriers 
are ineffective and suggests that the quantum dissociation mechanism is somehow bypassing the KAM barriers. 

A clue to the surprising quantum results comes from comparing the phase spaces in Fig. |ll| and Fig. |12| which show 
the uncontrolled and controlled cases respectively. Although, the KAM barriers seem to have reduced the extent of 
stochasticity, the 2:1 resonance is intact and appears to occupy slightly larger area in the phase space. Thus, this 
certainly indicates that a significant amount of the quantum dissociation is occuring due to the RAT mechanism 
involving the three Morse states n — 10, 12, and 14 as discussed before. In particular, the n = 12 state must still be 
actively providing a route to couple the initial state to the chaotic region via RAT. How can one test the veracity of 
such an explanation? One way is to scale the Planck constant down from h = 1 and monitor the quantum dissociation 
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FIG. 12. Phase space for the driven Morse system in the presence of the simpUfied classical control fields (cf. Eq. 

A2 ~ 0.01) designed to recreate specific KAM barriers (a) ujf /^r = 1 + 7^^ and (b) ujp/^r = \/3. Note that in both cases 

the dominant ^3,-2 Fourier amplitude of the leading order control field of Eq. |43] have been utilized and the desired KAM 

barriers are clearly seen (thick gray). The effect of the barriers seen in (a) and (b) on the classical (gray circles) and quantum 

(gray squares) dissociation probabilities are shown in (c) and (d) respectively. The uncontrolled results are indicated by the 

corresponding black symbols. Interstingly, the classical dissociation is suppressed but the quantum dissociation is slightly 

enhanced. 



process. Reduced Ti implies that the 2:1 island can support several states as well as the fact that other higher order 
resonances now become relevant to the RAT mechanism. Such a study is not presented here but one would anticipate 
that the dissociation mechanism can be understood based on the theory of RAT that already exists (See contributions 
by Schlagheck et al. and Backer et al., for example). Another way is to directly interfere locally with the resonance 
and see if the quantum dissociation is actually suppressed. Locally interfering with a specific phase space structure, 
keeping the gross features unchanged, is not necessarily a straightforward approach. However, the tools in the previous 
section allow for such local interference and I present the results below. 



The control field in Eq. 43 corresponding to the case of Fig. 12 'b) i.e., recreating the ujp/ilr — \/3 KAM barrier. 



turns out to be well approximated by 



Hc{J,e,t)KH{J,e,t)+ Y^ Fn^^2Cos{n9-2ujFt), 



(49) 



n=3,4 



The above Hamiltonian comes about due to the fact [45j that two Fourier modes i^3.-2 and -F4.-2 are significant in this 
case. Note that the specific KAM barrier of interest has a frequency between that of the 1+7^^ cantorus (around 



which the initial state is localized, cf. Fig. 11) and the 2:1 nonlinear resonance. From a perturbative viewpoint. 



creation of this KAM barrier using Eq. |49]is not expected to be easy due to the proximity to the 2:1 resonance and 
the fact that the F4__2 Fourier component is nothing but the 2:1 resonance. Nevertheless, Fig. 13 ^a) shows that the 
specific KAM barrier is restored and, compared to Fig. [Tlj the controlled phase space does exhibit reduced amount of 
chaos. Consisitently, Fig. lIsFc) shows that the classical dissociation is suppressed by nearly a factor of two as in the 
case shown in Fig. 12 h) wherein only the F^_2 component of the control Hamiltonian was retained. As mentioned 



earlier, in order to calculate the quantum dissociation probability the control Hamiltonian in Eq. 49 needs to be 
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FIG. 13. Phase space for the driven Morse system in the presence of the classical control fields designed to recreate the 
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Ldp/^r = \/3 KAM barrier. In (a) two dominant Fourier modes -F'3,_2 and -F4__2 of the leading order control term in Eq. 
retained. In (b) the simplified control field as in Eq. 47 is used with an effective field strength estimated using -f3,_2 and F4,_2. 
Notice that the desired KAM barrier is recreated in (a) but not in (b) (thick gray hue indicates the expected location). In (c) 
the classical dissociation probabilities of the uncontrolled (black circles), control using field (a) (open circles), and control using 
field (b) (gray circles) are shown. The quantum results are shown in (d) for the uncontrolled (black squares) and control using 
field (b) (gray squares) are shown. See text for discussion. 



mapped into a from as in Eq. |47| Since two Fourier modes need to be taken into account, the effective control field 
strength is given by 



F- 



3,-2 



F. 



4,-2 



V^{Jr) 14(Jr)' 



(50) 



Such a procedure yields A2 ~ —0.016 and thus the control field, still less intense than the primary field, comes with a 
relative phase of zero. Interestingly, as shown in Fig. 13 h), the resulting simplified control Hamiltonian fails to create 
the desired barrier. 

pre)). 



Moreover, the phase space also exhibits increased stochasticity and as a consequence the classical 
dissociation is enhanced (cf. 



Fig. 



However, Fig. 13 b) reveals an interesting feature - the 2:1 resonance is 



severly perturbed. This perturbation is a consequence of including the Fi-2 Fourier component into the effective 
control Hamiltonian. The key result, however, is shown in Fig. \l'^ d) where one observes that the quantum dissociation 
probability is reduced significantly. It seems like the quantum dynamics feels the barrier when there is none! The 
surprising and counterintutive results summarized in Fig. |12| and Fig. [13] can be rationalized by a single phenomenon 
- dynamical (resonance assisted in this case) tunneling. The main clue comes from the observation that quantum 
suppression happens as soon as the 2:1 resonance is perturbed. 

Although the results here are shown with a specific example, the phenomenon is general. Indeed computations (not 
published) for different sets of parameters have supported the viewpoint expressed above. Interestingly, the work of 
Huang et al. focused on classical suppression of ionization by creating local phase space barriers in case of the driven 
one-dimensional hydrogen atom^O]. Around the same time Brodier et al. highlighted [82] the importance of the RAT 
mechanism in order to obtain accurate decay lifetimes of localized wavepackets in the same system. Work is currently 
underway to see if the conclusions made in this section hold in the driven atomic system as well i.e., whether the 
attempt to suppress the ionization by creating local phase space barriers is foiled by the phenomenon of RAT. 
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C. Summary and future outlook 

The two examples discussed in this work illustrate a key point - Dynamical tunneling plays a nontrivial role in the 
process of quantum control. In the first example of the driven quartic double well it is clear that bichromatic control 
fails in regimes where chaos-assisted tunneling is important. In the second example of a driven Morse oscillator it 
is apparent that efforts to control by building phase space KAM barriers fail when resonance-assisted tunneling is 
possible. In both instances the competition between classical and quantum mechanisms is brought to the forefront. 
Although the two examples shown here represent the failure of specific control schemes due to dynamical tunneling, 
one must not take this to be a general conclusion. It is quite possible that some other control schemes might owe 
their efficiency to the phenonmenon of dynamical tunneling itself. Further classical-quantum correspondence studies 
on control with more general driving fields are required in order to confirm (or refute) the conclusions presented in 
this chapter. 

At the same time the two examples presented here are certainly not the the last word; establishing the role of 
dynamical tunneling in quantum coherent control requires one to step into the murky world of three or more degrees 
of freedom systems|83^. I mention two model systems, currently being studied in our group, in order to stress upon 
some of the key issues that might crop up in such high dimensional systems. For example, two coupled Morse oscillators 
which are driven by a monochromatic field already presents a number of challenges both from the technical as well as 
conceptual viewpoints. The technical challenge arises due to the fact that dimensionality constraints do not allow one 
to visualize the global phase space structures as easily as done in this chapter. One approach is to use the method of 
local frequency analysis[51] to construct the Arnol'd web i.e., the network of nonlinear resonances that regulate the 
multidimensional phase space transport. In a previous work [85', involving a time- independent Hamiltonian system, 
the utility of such an approach and the validity of the RAT mechanism has been established. However, "lifting" 
quantum dynamics onto the Arnol'd web is an intriguing possibility which is still an open issue. On the conceptual 
side there are several issues with multidimensional systems. I mention a few of them here. Firstly, even at the classical 
level one has the possibility of transport like Arnol'd diffusion[5Tl [72] which is genuinely a three or more DoF effect 
and has no counterpart in systems with less than three DoFs. Note that Arnol'd diffusion is typically a very long time 
process and is notoriously difficult to observe in realistic physical systemsjHB]. Moreover, arguments can be made for 
the irrelevance of Arnol'd diffusion (or some similar process) in quantum systems due to the finiteness of the Planck 
constant. Secondly, an interesting competition occurs in systems such as the driven coupled Morse oscillators. Even 
in the absence of the field the dynamics is nonintegrable and one can be in a regime where the modes are exchanging 
energy but none of the modes gain enough energy to dissociate. On the other hand, in the absence of mode-mode 
coupling, a weak enough field can excite the system without leading to dissociation. However, in the presence of 
such a weak field and the mode-mode coupling one can have significant dissociation of a specific vibrational mode. 
Clearly, there is nontrivial competition between transport due to mode-mode resonances, field-mode resonances and 
the chaotic regions[lSj. Selective control of such driven coupled systems is an active research area[HZ] today and the 
lessons learnt from the two examples suggest that dynamical tunneling in one form or another can play a central role. 
In the context of studying the potential competition between Arnol'd diffusion and dynamical tunneling, I should 
mention the driven coupled quartic oscillator system with the Hamiltonian: 

H{x,y,p^,Py;t) = -{pI + pI) + -{x"^ + y^) - ^ix'^y^ - xfoicosflit + coan2t). (51) 

In the absence of the field fo — the Hamiltonian reduces to the case originally studied by Tomsovic, Bohigas and 
UUmo wherein the existence of chaos-assisted tunneling was established in exquisite detail[TSl[T^. At the same time, 
with foy^O the system represents one of the few examples for which the phenomenon of Arnol'd diffusion has been 
investigated over a number of years. Nearly a decade ago, Demikhovskii, Izrailev, and Malyshev studied[88; a variant 
of the above Hamiltonian to uncover the fingerprint of Arnol'd diffusion on the quantum eigenstates and dynamics. 
An interesting question, amongst many others, is this: Will the fluctuations in the chaos-assisted tunneling splittings 
for /o = 0, observed for varying Ti, survive in the presence of the field? 

A related topic which I have not touched upon in this article has to do with the control of IVR using weak external 
fields. Based on the insights gained from studies done until now on field-free IVR (see also the contribution by Leitner 
in this volume), it is natural to expect that dynamical tunneling could play spoilsport for certain class of initial states 
that are prepared experimentally. The issue, however, is far more subtle and more studies in this direction might shed 
light on the mechanism by which quantum optimal control methods work. For instance, recent proposals on quantum 
control by Takami and Fujisaki[5n] and the so called "quantized Ulam conjecture" by Gruebele and Wolvnes[5D] take 
advantage (implicitly) of the system having a completely chaotic phase space. Dynamical tunneling is not an issue 
in such cases. However, in more generic instances of systems with mixed regular-chaotic phase space, a quantitative 
and qualitative understanding of dynamical tunneling becomes imperative. Given the level of detail at which one is 
now capable of studying the clasically forbidden processes, as reflected by the varied contributions in this volume, I 
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expect exciting progress in this direction. 
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